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Determining the structure of spectral densities is important for understanding the behaviour 
of any quantum field theory (QFT). However, the exact calculation of these quantities of¬ 
ten requires a full non-perturbative description of the theory, which for physically realistic 
theories such as quantum chromodynamics (QCD) is currently unknown. Nevertheless, it is 
possible to infer indirect information about these quantities. In this paper we demonstrate 
an approach for constraining the form of spectral densities associated with QFT propagators, 
which involves matching the short distance expansion of the spectral representation with the 
operator product expansion (OPE) of the propagators. As an application of this procedure 
we analyse the scalar propagator in (/^-theory and the quark propagator in QCD, and show 
that constraints are obtained on the spectral densities and the OPE condensates. In par¬ 
ticular, it is demonstrated that the perturbative and non-perturbative contributions to the 
quark condensate in QCD can be decomposed, and that the non-perturbative contributions 
are related to the structure of the continuum component of the scalar spectral density. 
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I. INTRODUCTION 


Spectral representations of matrix elements were first investigated by Kallen [l| and Lehmann j2|, 
and then later by and among others. An important consequence of these investigations was 
the discovery of the Kallen-Lehmann representation of the two-point function. For an arbitrary 
quantum field this representation relates the two-point function of the field (T{'I'(x)'I'(0)}) to 
an integral convolution between the free field propagator and some spectral density p. The integral 
representation enables one to determine interesting information about the analytic structure of 
correlation functions, and also has many important applications including the establishment of 
Goldstone’s theorem for relativistic local fields [^. Another important result in quantum field 
theory (QFT) is the operator product expansion (OPE). This expansion was first proposed by 
Wilson to describe the behaviour of products (or time-ordered products) of fields in the limit of 
coinciding space-time arguments. Given the renormalised fields A{x) and B{y), the OPE has the 
form: 


A{x)B{y)^Y.^,{x-y)Oi{y) (I.l) 

i=\ 


where {Oi(y)}(L^ is a finite set of renormalised fields, Ci are (possibly singular) coefficient functions, 
and ~ is understood to imply that an insertion of A{x)B{y) — Cli{x — y)Oi{y) into any Green’s 
function will vanish in the (weak) limit x —>• y. An important feature of the OPE is that both Oi{y) 
and the coefficients Ci depend on an auxiliary parameter p called the renormalisation scale. For the 
purpose of the discussions in this paper we are interested in the structure of two-point functions of 
certain fields T. By using the general form of the OPE outlined in Eq. m, these Green’s functions 
can be shown to have the following behaviour in the limit x ^ 0: 

(r{T(x)T(0)}) ~ Ciixmm (1.2) 


where (•) signifies the vacuum expectation value. The conceptual idea of the OPE is that the series 
provides an asymptotic decomposition of short and long distance degrees of freedom, which in 
the case of asymptotically free theories such as quantum chromodynamics (QGD) are partitioned 
between the Wilson coefficients C'i(x) and vacuum condensates (Cb(0)) respectively. For general 
theories though, this decomposition is not necessarily so clear-cut [7|. Nevertheless, the OPE has 
many important applications such as in the construction of factorisation theorems and the 
calculation of conformal field theories [^, as well as more applied uses like in the determination of 


QGD observables such as Rhad 


The spectral representation and the OPE are important results which have led to both successful 
experimental predictions and important theoretical developments. In particular, over the last 
few decades the determination of the perturbative and non-perturbative structure of QGD has 
significantly progressed due to the application of these results. The method that perhaps best 
epitomises the successful use of both the spectral representation and the OPE is the Shifman- 
Vainshtein-Zakharov (SVZ) sum rules ll|. By exploiting the analytic structure of certain correla¬ 
tion functions, this approach introduces a parametrised ansatz for the spectral density p and uses 
this to determine mesonic and hadronic parameters in terms of QGD vacuum condensates such as 
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(tpip) and Given lattice QCD estimates of these condensates, this then allows one to 

make a prediction for these parameters. The key point here is that it is not possible to exactly 
calcnlate the spectral density associated with a correlation function, the reason being that the 
complete analytic structure of QCD remains unknown. Instead, one has to constrain the form of 
p indirectly. Another example of a method which constrains the form of spectral densities is the 
so-called Weinberg sum rules [l^. These constraints are derived by performing a short distance 
expansion of the spectral representation of a correlation function, and inferring that certain linear 
combinations of the spectral densities must vanish if the correlation function in question has a 
specific singular behaviour. 

It is clear that constraining the form of the spectral density is very important if one wants to 
improve understanding of QCD, as well as other QFTs. In the literature this problem has been 
pursued in a variety of different ways, the SVZ and the Weinberg sum rules being two of the more 
developed methods. An interesting approach adopted by [I^ is to generalise the Weinberg sum 
rules by comparing the short distance spectral representation expansion of a correlator with its 
OPE. Based on which singular terms appear in the OPE, one can then conclude whether certain 
linear combinations of the spectral density vanish or not. In a similar manner, the authors in 
compare the expression generated by the large momentum propagator expansion in (/)'^-theory, with 
the leading singular terms in the OPE, but in this case with the intention of demonstrating the 
validity of the OPE itself. The success of this comparison approach between the short distance 
expanded spectral representation and the OPE, suggests that there may well be more information 
to be gained by performing a full expansion of both expressions, and then matching the resulting 
terms order by order in x. 

The remainder of this paper is structured as follows: in Sec. Ull we perform the short distance 
matching procedure for the scalar propagator in ^^-theory; in Sec. Illll we apply the same procedure 
to the quark propagator in QCD; and finally in Sec. |IV]we discuss the relevance of our results and 
the scope for further applications. 


II. SHORT DISTANCE MATCHING IN c^i^.THEORY 

In this section the short distance matching procedure outlined at the end of Sec. [I] will be applied 
to the propagator {T{(j){x)(j){0)}) in cj)'^ scalar field theory. Given the assumption of some standard 
QET axiomqj, this propagator has the following spectral representation: 

poo 

{T{(j){x)4>{0)})= ds p{s) iApix] s) (II.I) 

Jo 

where iAp{x] s) is the free bosonic Eeynman propagator, and p(s) is the spectral density. As with 
any QET, renormalisation of the fields is required in order to remove the divergences which arise 
as a result of the product of fields being ill defined at coincident space-time points. Once this 
procedure has been performed, the propagator instead satisfies the following renormalised spectral 


^ See [H-Q for more details. 
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representation 



{T{4)r{x)4ir{0)}) = / ds p{s,n,g) iAFix]s) 


(II.2) 


where (j)R is the renormalisation of the bare field cf), and the spectral density p is now also dependent 
on the renormalisation scale p and coupling g. If one now assumes x to be space-like < 0), 

the Lorentz invariance of the propagator enables one (for simplicity) to set 0- Under these 

conditions the free boson propagator iAj7’(x;s) has the following exact form 15l |: 


iApixo = 0, x; s) = 


dvr^lxl 


■Ki{y/s\x\) 


(II.3) 


where Ki is a modified Bessel function of the second kind. Under the assumption that the small-|x| 
behaviour of the integral in Eq. (jll.2p can be approximated by expanding the integrand around the 
point |x| =0, the propagator in this approximation is given by 


{T{4>r{x)(Pr{0)}) 


ds p{s,p,g) 


IGtt^ 


27-1 + 2 ln( Y ) + 2 ln(|x|) 


+ 


+ 0(|x|2) 


dvr^lxp 

Moreover, the renormalised propagator also has the following OPE [l^ : 

{T{(j)R{x)4>R{0)}) ~ Ci{x,p,m,g) + C^ 2 {x, p,m, g){(f)%{0)) + 


(11.4) 


(II.5) 


where m is the renormalised mass parameter, = [4’‘^]r is the renormalisation of the bare held 
0^, and • • • represents other possible non-singular terms. Under the assumption that x is space-like 
(with xq = 0), this asymptotic expansion is valid in the limit |x| —?■ 0. The Wilson coefficients 
Cl and C ^2 can be calculated perturbatively, and it turns out that to lowest order in perturbation 
theory they have the following form 17j: 


MpM^)+0{g^) 


1 Tfl 

C^ 2 {x,p,m,g) = 1 + (+Vn +C’(5'^) 


(11.6) 

(11.7) 


Inserting these expressions into Eq. (111.5^ then gives 


{T{(l)R{x)(j)R{0)}) 


+ 


m 


dvr^lxp IGtt^ 


In (/iVP) + In ('/>l(0) + (0|(O)) + 0{g‘^) 

{11.8) 


Since Eqs. (III.4P and (jll.sp correspond to equivalent descriptions of the propagator in the small-|x| 
limit, and the spectral density p is not x dependent, one can equate these two equations and match 
the coefficients of the various |x|-dependent terms. In doing so, one obtains the following relations 
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between the OPE coefficients and certain moments of the spectral density: 

ds pis,fi,g) = l (II.9) 

1*00 

O (|x|°) : / ds sp{s,p,g) 

Jo 

+ 5ln(^)(</.|(0)) + O(ff2) (II.IO) 

0 ( 111 ( 1 x 1 )):/" ds sp{s,p,g) = m^+ +0{g‘^) (11.11) 

One thing to notice here is that the relation in Eq. pL9p is exact since there are no other O 
terms in Eq. pi.4p . and it is not possible to generate another term with this dependence in Eq. pi.Sp 
no matter what perturbative order C ^2 and Cj are expanded to. Equations (lll.iop - dll.lip on the 
other hand are only perturbatively valid to 0 ( 5 ^), since expanding 0^2 and On to higher orders 
may generate additional constant or 0(ln(|x|)) terms. By inspecting Eq. pi.lip , it is clear that 
this is satisfied if the spectral density is given by 

p{s,p,g) = S{s - m^) + |(</>h(0))-4(s) + 0{g^) (11.12) 

where ^(s) satisfies the normalisation condition 

ds sA(s) = 1 (11.13) 

and also implicitly depends on the renormalisation scale p. Erom Eq. pi.l2p one can see that the 
spectral density has several interesting features: there is a Dirac delta term which corresponds 
to the existence of a state with mass m in the theory; the second term has the structure of a 
continuum component since it contains an explicit factor of the coupling constant g, and is hence 
a by-product of interactions in the theory; and also the second term is premultiplied by the con¬ 
densate (</>^(0)), which suggests that the contribution of this continuum component to the spectral 
density is moderated by the magnitude of the scalar condensate. 




27 - 1 -h 2 In 


= 2m^ In (p) -|- 167r^((/>|.(0)) 


Inserting this expression for the spectral density into Eq. pi.101) . and ignoring terms of 0{g‘^) and 
above, one obtains the relation 


|(7r(0)) [27-1.1(4) 


l+X\ + m^ 


27 


1-F2 In 



2w? In {p) -|- 167r^((/>|j(0)) 

+ g\n{p){(fim+0{g^) (11.14) 


where X has the form 


Upon rearrangement this gives 


X = 


poo 

/ ds sln(s)74(s) 
Jo 


(11.15) 




g 


Ibvr^ 

2 


X' 


-1 


C-^lnfii 

Svr^ \m 


+ 0{g^ 


Svr^ \m/ IGvr^ 
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where X' and C are 


In (/t) - ^ [27 - In (4) - 1 -h X] 

(11.17) 

9 

TTl 

(11.18) 


A significant feature of the expression for the condensate in Eq. pi.161) is that it explicitly depends 
on X, an integral involving the a priori unknown continuum contribution to the spectral density 
A(s). Because this condensate does not receive any non-perturbative contributions 0], it must 
have exactly the same form as the purely perturbative expression for (p^(0)) computed using the 
renormalisation of the operator ^^(0). Therefore, by equating these expressions one can obtain 
information about the continuum component A(s). 


In general, a renormalised operator satisfies the following renormalisation group equation 
(RGE) 0: 


3 


(11.19) 


where is the anomalous dimension matrix and {Ojj} is a finite closed basis of renormalised 
operators with dimension < dim(Cl)^). In i^'^-theory </>|j(0) mixes with the identity operator I, but 
not with 16]. The RGE for therefore has the form 


d 


— 7(^2^2())|j( 0) +7^2jl 


( 11 . 20 ) 


By definition, the vacuum expectation value of p^(0) with the perturbative (Eock space) vacuum 
state vanishes. However, the vacuum expectation value with the physical non-perturbative vacuum 
state does not necessarily vanish, and this is what (p^(0)) corresponds to both in the preceding 
and proceeding discussions in this section. After inserting both sides of Eq. (III.20p between the 
physical vacuum state, one obtains the following RGE for (i?i^(0)): 


where /3 is the beta function of the theory and 7 ^ is the anomalous mass dimensiorH. At one-loop 
order one has 0,0 


13 = 


lEL 

167r^ ' 


70^0^ 27m 


167r^ ’ 


7021 


m 


( 11 . 22 ) 


By choosing a mass-independent renormalised operator basis, in this case the anomalous 

dimensions can in general become mass dependent 1^, and this is in fact what happens for 7 (^ 2 ][. 


Here we use the opposite sign convention to [1^ for 7 ^ 
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Using the method of characteristics, the solution of Eq. (III.2ip is equivalent to the solution of the 
following set of ordinary differential equations; 


dln/i 
dt ~ ^ 

(11.23) 

d9 _ o_ 3^^ 
dt IGvr^ 

(11.24) 

dm gm 

dt = = 32vr2 

(11.25) 

(</’r(0)) - 7</2</2 (Pi?(0)) -1- 7</2i - (<(’ij(0)) g^2 

(11.26) 


However, in order to obtain unique solutions one must first specify a boundary condition for each 
equation. Since the variable t has no physical significance and only serves to parametrise the 
characteristic curves along which solutions are defined, one can choose all of the boundary data 
to be at t = 0. For Eq. pi.23p the general solution is given by In// = t + ci, so the integration 
constant has the form ci = ln//(t = 0). Letting ci = ln(/2), where fl is some physical scale, the 
condition t = 0 is equivalent to // = //, and so t = In With this choice of boundary condition 
the solutions of Eqs. pi.24l) - (III.25l) can be written 


9 = 9 



m = m 



IGvr^ ) 


where g = g{t = 0), m = m{t = 0), and the solution of (p|j(0)) has the form 


(4(0)) 


2m? 

9 



2m^ 

9 



+ (4(0)) 



(11.27) 


(11.28) 


with ((/>^(0)) = (0^(O))(t = 0). By inverting the expressions in Eq. (III.27h . one can rewrite the 
solution in Eq. (III.281) exclusively in terms of the parameters g, m and t = In (Doing so gives 


2?T7 ^ _ 

(4(0)) = — + (4(0)) 

9 


1 + 


3g 

IGvr^ 


ln(^ 

9 


2m^ 


1 + 


3g 

IGvr^ 


ln(^ 

9 


(11.29) 


Because this perturbative determination of (4(0)) is valid up to one-loop order, the solution is 
therefore equal to the following expansion of Eq. (III.291) up to 0{g)\ 


(«(0)> = («(0))-^ln(^)-^ 


In 


9 


gm 




In 


9 
9J \ 


+ 0{g^) (11.30) 


Finally, one can now compare this equation with the expression for (4(0)) [(Eq. pi.lGI) ] obtained 
via the spectral density matching conditions in Eqs. (III.9p - (III.lll) . One can clearly see that these 
expressions have a very similar form. In fact, using the solutions for g and m, one can rewrite 
Eq. (III.161) as follows: 


m 


(4(0)) = ^ ^ 




IGvr^ 


IGvr^ 


X'C + 


gm 

1287r4 


C In 


\n{^]X' + 0{g^ 


gm^ 

+ —7 In 


25677 ^ 


m 


(11.31) 
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where the constant C is defined as 


C = 


_ 9 

m 

IGvr^ 


[ 27 -In (4)-1] 


By demanding that X' satishes the following relation 


I' = ln(4) + 


21n(A)h-ln(2)] 

m/ ' [[2^-l„(4)-ll-21n(a] 


the expression for (i?i^( 0 )) becomes 






IGvr^ \mJ 1287r^ 1 \mJ. 
which has exactly the same form as Eq. (III.30p if one makes the identihcation 


(11.32) 


(11.33) 


(11.34) 


= {<Pr{0)) = C (11.35) 

So equating the short distance matched and RGE derived expressions for (<?i>^(0)) has introduced 
two new constraints: the functional form of the initial conditions in Eqs. (III.23P ~ (III.26P is hxed, 
and hence the form of (</>^(0)) is completely specified in terms of the free parameters m and g; and 
the condition in Eq. pi.33p implies that ^ 4 ( 5 ) must satisfy 



ds s ln(s)24(s) = ln(4m^) — 27 + 1 


41n(ah-ln(2)] 

[|27-l.i(4)-l]-21.i{i)] 


and therefore provides an additional constraint on the form of the spectral density p. 


(11.36) 


Although (/)'^-theory may well not be physically realistic due to its triviality 19|, [20|] , the discussion in 
this section demonstrates that the short distance matching procedure provides a way of determining 
new constraints and qualitative features of the theory, and in particular the spectral density, which 
contrasts with numerical-based approaches 


21 - 2 , 


Moreover, because this procedure is model 
independent, since it only relies on the existence of an OPE and a spectral representation, it can 
also equally be applied to physically realistic theories such as QCD, and this is what we pursue in 

Sec. in 


III. SHORT DISTANCE MATCHING IN QCD 


The short distance matching procedure that was performed for </>'^-theory is equally applicable 
to QCD, and in this section we focus in particular on analysing the fermionic quark propagator 
{T{ipix)'ip(0)}) in this way. For this propagator, the spectral density p can be decomposed in spinor 


space as 


m 


p{s) = ^ 5 ( 5)1 + PPs(s)75 + Pvishfi + Pp\/('S)757m + Pt (III.l) 


where the spinor indices are suppressed. It turns out that the tensor term in Eq. (jlll.lh does not 
contribute, and furthermore if one assumes the absence of parity violation, then pps = ppv = 0. 

























Combining these results, the quark propagator has the following renormalised spectral representa¬ 
tion: 

POO 

{T{'iI>r{x)^r{0)}) = / ds pvis,iJ,,g) iSFix-,s)-\-iAF{x-,s)[psis,iJ,,g) - y/7pv{s,n,g)] (III.2) 

Jo 

where iAF^x] s) and iSF{x] s) are the free bosonic and fermionic Feynman propagators respectively, 
= P^) '■= P^Pv{s), and ipR, 'ipR are the renormalised bare helds. Assuming x is space-like (and 
setting xo = 0), one can perform a small-|x| expansion in an analogous way to Sec. [Ill The space-like 
structure of the free bosonic propagator is given by Eq. ()II.3p . and for the free fermionic propagator 
it has the form 


iSF{xo = 0, x; s) = [{i^ + V7) iAF(x; s)] 

's [Ko{^/s'\x\) + K2{^/F\x\)] y/s 


= -^7 Xi 


+ 


87r^|x|' 




47r^|x| 


;^l(v^|x|) 


dvr^lx 


■Ki{y/s\x\) 


(III.3) 


Finally, inserting the explicit expressions for the free propagators into Eq. (IIII.2p . and expanding 
around the point |x| = 0, one obtains 


POO 

{T{iI!r{x)^r{ 0)})^ dspv{s,p,g) 

Jo 




27r^|x|^ Svr^lxp 


+ ps{s,p,g) 


IGtt^ 


2 In 


(^ 1 + 27-1 + 2 ln(|x|) 




47r^|x|- 


+ 0(|x|2) 


(III.4) 


In a similar manner to (/>^-theory, the quark propagatoiH has the following operator product expan- 

(III.5) 


Sion 


2^: 


{T{iPr{x)iPr{&)]) ~ Ci{x,p,m,g) + C:;p^{x, p,m, g){'tp-ip{0)) H- 

where Ci and satisfy the RGEs 

d d 

~ ^~dp^^ ~ (III.6) 

With this RGE convention for the Wilson coefficients, the anomalous dimensions at one-loop order 
are given by 130 


9 


3m'^ 


27r2 ’ 271^ I ^ 


9 


(III.7) 


® For simplicity we assume here that there is only one flavour of quark, with mass m. 

As in Sec. El we adopt a mass-independent renormalised operator basis here (like [ 2 rI]'I. which means that the 
anomalous dimensions can in general be mass dependent, unlike in [^. Nevertheless, the mass-dependent anoma¬ 
lous dimensions are related to the mass-independent ones by a multiplication of a certain power in the mass m (in 
this case 7 ^,^ j = these choices lead to the same RGE for tp'iplO). 
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Just like in Sec. mi one can solve these equations using the method of characteristics. In this case, 
solving these equations (to one-loop order) requires one to solve the following ordinary differential 
equations: 


ding ^ 
dt 

(III.8) 

dg _ o_ ^g^ 

dt lOvr^ 

(III.9) 

dm g‘^m 

dt = = 2vr2 

(III. 10) 

d r. _ 9^ ^ 

dt^ 27r2 w 

(III.ll) 

d / o2 X 

dt^ 27r2 1 +37r2j^^V> 

(III.12) 


With the initial conditions n{t = 0) = g{t = 0) = g, and m{t = 0) = m, one has 


2 -2 I 1 , 

0 =» 1^+8? 


m = m 1 -|- 


7ft 

8^ 


(III.13) 


where t = In (/r|x|), and the Wilson coefficients have the following form; 




- 1 - 


7gH 

8^ 


2.x I 


Ci = - 




27r2|x|^ 
4m^ 




m 




dvr^lxp 

,2 


iilm? Arf 
H- - -1- 

^ . n, ,o ^ 3^2 

2 . 


SvT^lxp 

1 + + 14t) (l - g) 


1 + 


^f 


1 - 


7gh 


- 1 ' 


IGvr^ \ Stt^ 
7gH 


- 1 ' 


.2.\ f 


1 - 


87r2 


(III. 14) 


1 - 


87r2 


2 + \ 7 


(III.15) 


Expanding these expressions to 0{g^) and inserting them into Eq. (IIII.5P gives 


(r{'0R(x)V’R(O)}) ~ - 






m 


27r2|x|^ dvr^lxp Svr^lxp 27r2 


iitm? 3m^ . / , ix g'^rn? , , , 

+ . oi lo + In (/r|x|) In (Ai|x|) 


+ [In - W(0)) + ^ In (/x|x|) (V’V'(O)) + 0{g^) 


(III.16) 


Since both the spectral densities are x independent, one can perform the same procedure as in 
Sec.lH and match the different |x|-dependent coefficients in this expression with the moments of 
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the spectral density in Eq. IIII.4t 

roo 

ds pvis,fi,g) = 1 


o{ 


Wr J 

o{ 


WrJ 

o( 


\xpj 

0(1 

x|0) 


ds spv{s, p, g) = m? 

0 

poo 

(|x|°) : / ds sps{s,p,g) 

Jo 


(III. 17) 

(111.18) 

(111.19) 


27 - 1 + 2 In ( 


= 24m^ In (p) + 


8g^m^ 


TT^ 


In (p) 


+ l^^[ln(^)]2_l67r^(V;V.(0)) 


vr^ 


O 


(ln(|x|)) : f 

Jo 


ds sps{s, p, g) = 12 m^ + 


Ag‘^m^ 


+ 


+ 8g‘^ In {p) (V'V’(O)) + 0(5 
125 r^m^ 


2(V-- 

S', 


(III.20) 


vr^ 


TT^ 


■ln{p)+4g^{i^i;{0))+O{g^) ( 111 . 21 ) 


From Eq. piI.21D one can see that this relation is satisfied if the scalar spectral density has the 
following form: 


Ps(s,p,g) = 


Ag"^ 12g^ 

12 H- 5 —I-^ In (/i) 


m6{s — m^) + 4g^(V^'i/’(0))-6('S) + 0 ( 5 ^) 


where B{s) satisfies the normalisation constraint 

ds sB{s) = 1 


/ 


(III.22) 


(III.23) 


and also implicitly depends on the renormalisation scale p. It is interesting to note here that ps 
has the same characteristics as the spectral density in Eq. pi.l2p : a Dirac delta term, and a 
continuum contribution B{s) which has an explicit coupling constant and condensate prefactor. 

Similarly to Sec. [Ill by substituting ps into Eq. pil.20p one can rearrange to obtain an explicit 
expression for the quark condensate: 


(V^V^(O)) = /C + ^ln 


+ 


2Tr‘^ 
8g^m^ 


fp\, g-m'^ ^ ( P\ , 9“ , 8g‘^m? ^ ( P \ 

Km) + ^ (m) + ^ (m) ^ 

+ 0{g^) 


[ln(^)] — 2 In (m) In (/r) 


+ /C 


^2 




where JC and JT' are given by 


J' = ln{p)--[2^-ln{4)-l+J] 
/C = -^[27-ln (4)-l] 




and 77 is dehned as 


J = 


f 


ds s In (s) B{s) 


(111.24) 

(111.25) 

(111.26) 

(111.27) 
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In just the same way, this condensate explicitly depends on the unknown continuum component 
of the spectral density B{s). However, unlike the scalar condensate in (/>^-theory, {’ip'ip(0)) contains 
both perturbative and non-perturbative contributions 2^, and as we demonstrate it turns out that 
the non-perturbative contributions arise due to B{s). To make this more precise, one must first 
calculate the perturbative contributions to ('^'(/^(O)) [denoted (V^V’(0))p] which originate from the 
renormalisation of '0V'- The RGE of (00(O))p is 


d 




(III.28) 


where 7^^^^ and 7^^]i are given in Eq. (jTTTTj). With the boundary condition (00(O))p(t = 0) = 
(00(O))p, the solution to Eq. pil.28p has the form 


(00(O))p = 


4m^ 


4m^ 

~w 


1 + 


3g" 


1 - 


7gH 

87r2 


- 1 ' 


1 - 


7gH 


87r2 
, 2 .\ -1 




1 - 


(00(O))p 

IgH 


1 - 


JgH 

Svr^ 


Stt^ 


(III.29) 


Because the anomalous dimensions used in Eq. (IIII.28h are only valid up to O(g^), the perturbative 
expansion of (00(O))p is also only valid up to this order. Performing this expansion gives 


(00(O))p = (00(O))p + 


3m^ 


27r2 


In ( ^ 1 + 


2 

g 




g 


+ |^ln(^) W(0))P + 


3g^m^ ^ 
27r4 


ln(^ 

g 


Inl^ 

g 


+ 0(g3 


(III.30) 


By also expanding Eq. (IIII.24I1 to 0{g^), and comparing this expression with Eq. pil.30p . the quark 
condensate can be decomposed as follows: 


(00(0)) = (00(O))p 


27r2 


C + ?01n 


27r2 


+ 


ii 

fn 

— m2 


In (2m) - 7 + -(1 - J") 


2 [In (m)]" - [27 - In (4) - 1] [3 hi (i) - I - 2 In (m)] - 2 


[[ 27 -ln( 4 )-l]- 21 n(a] 


(III.31) 


where g = m and (00(O))p has the form of Eq. pil.301) with 




3m 


(00(O))p =JC = [27 - In (4) - 1] 


(III.32) 


Since the first term is purely perturbative, it must be the case that the second term parametrises the 
non-perturbative contributions to the quark condensate, and in particular, the integral involving 
B{s). This explicit decomposition of the quark condensate into perturbative and non-perturbative 
contributions has not to our knowledge been established before in the literature, and instead has 
simply been assumed 27l|. Moreover, the direct connection between the non-perturbative contribu¬ 
tions and the continuum component of the scalar spectral density B(s) has not been made before. 
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This has interesting applications because it means that if one can estimate the form of B(s) from 
the integral constraints in Eqs. (IIII.18D . (IIII.20p and (IIII.23D . one can use Eq. piI.31D to directly 
estimate the non-perturbative component of . 


The analysis in this section has demonstrated that by equating the short distance expansion of 
the spectral representation of the quark propagator in QCD with its OPE, one can obtain novel 
information. A nice feature of this method, by contrast to more numerical-based approaches 28l43fll| . 
is that it requires practically no theoretical input other than the form of the Wilson coefficients, 
and yet from this one is able to derive the qualitative structure of the scalar spectral density 
psi impose integral constraints on both ps and pv, and explicitly decompose the perturbative and 
non-perturbative contributions to {ip'ifjiO)). Moreover, unlike techniques such as the SVZ sum rules, 
phenomenological approximations such as quark-hadron duality 3l[ are not assumed, which makes 
this approach process independent and therefore applicable to arbitrary correlators. In principle, 
this approach could also provide useful input for the SVZ sum rules. A key feature of these sum rules 
is the requirement to introduce a parametrised form of a spectral density [s^, and so information 
obtained about the structure of this spectral density from the short distance matching procedure 
could be used to provide additional constraints on the corresponding parameters. 


IV. CONCLUSIONS 


Spectral densities play a central role in determining the dynamics of a QFT, and yet in many 
instances it is not possible to calculate these objects exactly. This obstruction arises because the 
non-perturbative structure of these theories is not well understood. Nevertheless, one can infer 
information about the form of spectral densities by applying general QFT techniques. In partic¬ 
ular, in this paper we have demonstrated that by matching the short distance expansion of the 
spectral representation of the scalar propagator in (/>^-theory and the quark propagator in QCD 
with their respective OPEs, constraints on both the spectral densities and the OPE condensates 
arise. On a qualitative level these constraints are interesting because they provide new information 
about the form of the spectral densities, and specifically the structure of the continuum contribu¬ 
tion. In the case of QCD, this information can then be used to explicitly decompose the quark 
condensate into perturbative and non-perturbative contributions, and it turns out that 

the non-perturbative contributions are related to the structure of the continuum component of 
the scalar spectral density. More directly, these constraints may also provide useful information 
for procedures such as the SVZ sum rules which rely on constructing a parametrised form of the 
spectral density of certain correlation functions. A nice feature of this short distance matching 
approach is that it is completely model independent - it only relies on the existence of an OPE and 
a spectral representation. So in principle the analysis applied to the scalar and quark propagators 
in this paper can equally be applied to other interesting correlators such as the gluon propagator, 
the vector current correlator (T {J^(x)Jiy{0)}), or other more general matrix elements, and this 
could potentially provide some interesting new insights. 


13 














ACKNOWLEDGMENTS 


I thank Thomas Gehrmann for useful discussions and input. This work was supported by the Swiss 
National Science Foundation (SNF) under contract CRSII2_141847. 


[1] G. Kallen, Helv. Phys. Acta 25, 417 (1952). 

[2] H. Lehmann, Nuovo Cimento 11, 342 (1954). 

[3] R. Jost and H. Lehmann, Nuovo Cimento 5, 1598 (1957). 

[4] F. Dyson, Phys. Rev. 110, 1460 (1958). 

[5] F. Strocchi, Symmetry Breaking (Springer-Verlag, 2008). 

[6] K. Wilson, Phys. Rev. 179, 1499 (1969). 

[7] V. A. Novikov, M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B249, 445 (1985). 

[8] J. Collins, Foundations of Perturbative QCD (Cambridge University Press, Cambridge, 2011). 

[9] P. Deligne, P. Etingof, D. S. Freed, L. C. Jeffrey, D. Kazhdan, D. A. Morrison, and E. Witten, Quantum 
Fields and Strings: A Course for Mathematicians^ Vol. 1 (American Mathematical Society, Providence, 
1999). 

[10] K. G. Chetyrkin and J. H. Kiihn, Nucl. Phys. B432, 337 (1994). 

[11] M. Shifman, A. Vainshtein, and V. Zakharov, Nucl. Phys. B147, 385 (1979). 

[12] S. Weinberg, Phys. Rev. Lett. 18, 507 (1967). 

[13] C. Bernard, A. Duncan, J. LoSecco, and S. Weinberg, Phys. Rev. D 12, 792 (1975). 

[14] S. Schweber, An Introduction to Relativistic Quantum Field Theory (Row, Peterson and Company, 
Evanston, 1961). 

[15] K. Huang, Quantum Field Theory: From Operators to Path Integrals (John Wiley & Sons, New York, 

2010 ). 

[16] J. C. Collins, Renormalization (Cambridge University Press, Cambridge, 1984). 

[17] H. Osborn, Encyclopedia of Mathematical Physics: Operator Product Expansion in Quantum Field 
Theory, Vol. 3 (Elsevier Ltd., New York, 2006) p. 616. 

[18] H. Kleinert and V. Schulte-Frohlinde, Critical Properties of (j)^- Theories (World Scientific, Singapore, 

2001 ). 

[19] J. Glimm and A. Jaffe, Quantum Physics: A Functional Integral Point of View (Springer-Verlag, 1987). 

[20] F. Strocchi, An Introduction to Non-Perturbative Foundations of Quantum Field Theory (Oxford Uni¬ 
versity Press, Oxford, 2013). 

[21] M. Jarrell and J. E. Gubernatis, Phys. Rep. 269, 133 (1996). 

[22] Y. Burnier and A. Rothkopf, Phys. Rev. Lett. Ill, 182003 (2013). 

[23] D. Dudal, O. Oliveira, and P. J. Silva, Phys. Rev. D 89, 014010 (2014). 

[24] C. Itzykson and J. B. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1980). 

[25] F. J. Yndurain, The Theory of Quark and Gluon Interactions (Springer-Verlag, Berlin, 1999). 

[26] K. G. Chetyrkin and A. Maier, J. High Energy Phys. 01, 92 (2010). 

[27] V. P. Spiridonov and K. G. Chetyrkin, Sov. J. Nucl. Phys. 47, 522 (1988). 

[28] R. Alkofer, W. Detmold, C. S. Fischer, and P. Maris, Phys. Rev. D 70, 014014 (2004). 

[29] F. Karsch and M. Kitazawa, Phys. Rev. D 80, 056001 (2009). 

[30] S. X. Qin, L. Chang, Y. X. Liu, and C. D. Roberts, Phys. Rev. D 84, 014017 (2011). 

[31] M. Shifman, At The Frontier of Particle Physics (World Scientific, Singapore, 2001) Chap. 33, p. 1447. 

[32] P. Colangelo and A. Khodjamirian, At The Frontier of Particle Physics (World Scientific, Singapore, 
2001) Chap. 34, p. 1495. 


14 



